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RADIATION FROM A SLOT ANTENNA IN A GROUND PLANE 


COVERED BY A WARM PLASMA LAYER* 
by Norman C. Wenger 
Lewis Research Center 

SUMMARY 

The electromagnetic properties of a model of a reentry plasma layer are determined 
by using the linearized hydrodynamic equations for a single fluid warm plasma medium. 
The model consists of a perfectly conducting plane of infinite extent covered by a warm 
plasma layer of uniform density. The field excited by a slot antenna in the conducting 
plane is computed by using Fourier transform techniques. Expressions are obtained for 
the radiation field and for the waves guided along the plasma layer. Numerical results 
are presented that show the dependence of the energy transport by the field on frequency 
and electron temperature. In addition, the radiation pattern of the far-zone field is com- 
puted for various frequencies and electron temperatures. 


INTRODUCTION 

A considerable amount of research has been done in recent years on the electromag- 
netic aspects of space flight. The particular subject that has perhaps received the great- 
est attention is the effect on communication systems of the ionized gas layer that envelops 
a reentry vehicle. In a typical reentry situation, the ionized gas layer attenuates radio 
signals to such an extent as to cause a communications blackout. It is, therefore, of 
prime importance to understand reentry phenomena as completely as possible in order 
to develop methods to alleviate the blackout problem. 

The proposed solutions of the blackout problem can generally be grouped into two 
categories. In the first group are those methods that attempt to alter the ionized 
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Portions of this report were included in a thesis submitted to Case Institute of 
Technology, Cleveland, Ohio, May 1965 in partial fulfillment of the requirements for the 
degree of Doctor of Philosophy. 



gas layer by injecting a chemical into it in order to produce electron attachment or to 
hasten the recombination processes (ref. 1). The second group contains those methods 
whereby a magnetic field is impressed upon the ionized layer in order to make it more 
transparent to electromagnetic radiation (ref 2). Before the effectiveness of any of 
these methods can be determined, it is necessary to have a valid model of the ionized 
layer. 

The majority of the theoretical analyses on the effect of an ionized gas on electromag- 
netic radiation use the cold plasma model. A common configuration employed in the 
treatment of the reentry problem consists of a perfectly conducting plane covered by a 
cold plasma layer of uniform density. The source for the electromagnetic radiation is 
usually a slot antenna in the conducting plane. An extensive account of the major phenom- 
ena observed for this model appears in references 3 to 5. The results of these investiga- 
tions indicate that the cold plasma layer has a very small effect on the electromagnetic 
radiation for frequencies sufficiently above the plasma frequency. At frequencies near 
or below the plasma frequency, severe distortion or attenuation of the radiation field is 
present. In addition, there may be waves trapped by the plasma layer. Some of these 
trapped waves exhibit very unusual properties such as being backward waves. A knowl- 
edge of the properties of these waves is extremely important in determining the efficiency 
of an antenna system since these trapped waves may transport a considerable amount of 
energy. 

One of the main criticisms that can be made against the use of the cold plasma model 
in connection with reentry studies is that the very mechanism that produces the plasma, 
namely, high temperature, has been neglected. Some work has been done that uses the 
warm plasma model. Stewart and Caron (ref. 6) have computed the radiation from a line 
source in a ground plane covered by a warm plasma layer. The results of their analysis 
are, however, limited by an approximation that requires the source frequency to be near 
the plasma frequency. 

This investigation presents an exact analysis of a model of a warm reentry plasma 
layer to determine the effect of the plasma temperature. The total field excited by a slot 
antenna in a conducting plane covered by a warm plasma layer is computed by using 
Fourier transform techniques. Expressions are obtained for the far-zone radiation field 
and for the waves guided by the plasma layer. The energy transported by the radiation 
field and by each of the guided waves is computed. The spectrum of waves associated 
with a plasma layer is shown to be significantly changed by the nonzero plasma tempera- 
ture. The effect of the temperature is to generate a completely new spectrum of guided 
waves that are electroacoustic in nature in addition to perturbing the waves associated 
with the cold plasma layer. 

This research was partly supported by the Air Force Cambridge Research Laborato- 
ries, Office of Aerospace Research, under contract AF 19(628)-1699. 
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SYMBOLS 


a 


a , a , a 
x’ y’ z 


b 

C 


c 

D(0) 

d 

E 

e 

F 

if 

J 


h 

h 


y 


Im 


I 


m 


acoustic velocity in electron gas 
unit vectors in x-, y-, and z- directions 

plasma layer thickness 

inversion contour (fig. 2) 

velocity of light 

see equation (35g) 

width of slot antenna 

total electric field intensity 

perturbation electric field intensity 

long-range force acting on electrons 

total magnetic field intensity 

Fourier transform of h 

perturbation magnetic field intensity 
y- component of h 

imaginary part of 
magnetic current 


total electric current density 
total magnetic current density 

jjH perturbation magnetic current density 


k Boltzmann constant 

k 0 wave number for magnetic field in plasma 

k Q free-space wave number 

k wave number for pressure field in plasma 

Xr 

M(x, / 3 ) see equation (35f) 

m mass of electron 
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N 


total electron number density 


N(x, /3) see equation (3 5e) 
n perturbation electron number density 

n Q equilibrium electron number density 

P total electron pressure 

cP Fourier transform of p 

p perturbation electron pressure 

p Q equilibrium electron pressure 

q electronic charge 

Re real part of 

r radial coordinate (fig. 1) 

r general position vector 

S energy flux vector 

T g electron temperature 

t time 

u(z) unit step function 

V total electron fluid velocity 

v perturbation electron fluid velocity 

x, y, z rectangular coordinates (fig. 1) 

a f ^ 

(3 Fourier transform variable 

n** 1 root of D(/3) = 0 

r h transverse wave number for magnetic field in plasma 

r p transverse wave number for pressure field in plasma 

y ratio of specific heats with constant pressure to constant volume of electron gas 

6(x) delta function 

e Q electric permittivity of free space 
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e relative dielectric constant of plasma 

ir 

rj imaginary part of v 

9 polar coordinate 

9 critical angle associated with magnetic field 

9 q critical angle associated with pressure field 

Xg wavelength of guided waves in plasma 
X Q free- space wavelength 

X free-space wavelength at plasma frequency 

Jr 

/i Q magnetic permeability of free space 

v complex variable 

p modulus of v - (7t/2 - 9) 

a real part of v 

<p phase of v - (n/2 - 9) 

u> source frequency 

plasma frequency 

r 

* complex conjugate 


BASIC EQUATIONS FOR A WARM PLASMA 

A plasma can be considered as being comprised of several miscible fluids. These 
fluids will consist, in general, of an electron fluid, an ion fluid for each species of ion, 
and a neutral particle fluid for each species of neutral particle. The dynamics of these 
fluids will be described by the basic hydrodynamic equations. An excellent discussion 
on the validity of the hydrodynamic treatment has been given by Oster (ref. 7). 

In order to simplify the analysis of a warm plasma, it will be assumed that the time 
variation of all nonstationary quantities is sufficiently rapid so that only the electron 
motion need be considered. This assumption allows the plasma to be described by the 
single electron fluid model. For this case, the transport equations of interest are 
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( 1 ) 


+ V • (NV) = 0 
3t 


mN — + mN( V • v) V = NF - VP 

at 


( 2 ) 


PN - ^ = constant 


(3) 


where 

N(r, t) electron number density 

V( r, t) electron fluid velocity 

P(r, t) electron pressure 

F( r, t) long-range force per electron 

y ratio of specific heat with constant pressure to specific heat with constant vol- 

ume of electron gas 

m mass of electron 

Equations (1) to (3) are a statement of the conservation of matter, conservation of mo- 
mentum, and conservation of energy of the electrons, respectively. Dissipative effects 
such as collisions between the electrons and ions have been neglected to simplify the anal- 
yis. 

Maxwell’s equations must be employed to complete the description of the dynamics of 
a warm plasma: 


V X H 


J + £ — 

° at 


(4) 


VXE=-J m 



(5) 


The new symbols used are 

J( r, t) electric current density 
J m ( r, t) magnetic current density 

H( r, t) magnetic field intensity 
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E( r, t) electric field intensity 
e Q electric permittivity of free space 

magnetic permeability of free space 

The electromagnetic equations have been formulated so that the electric current density 
J includes all conduction, convection, and polarization currents in the plasma. A mag- 
netic current density J m has been included in equation (5) since it provides an easy way 
to introduce the effect of the slot antenna later. 

The hydrodynamic equations and the electromagnetic equations are coupled by the 
force F and the current density J. If gravitational fields are neglected, F can be given 
by the Lorentz force equation 


F = q( E + V X jLt Q H) 

where q is the charge on the electron. The current density J is given by 


( 6 ) 


J = qNV (7) 

for the case where the conduction and polarization currents do not exist. This situation 
will prevail if there are no electric current sources within the plasma and if the polariza- 
tion field of the neutral particles and ions is negligible compared with the applied field. 

Equations (1) to (7) represent 17 scalar equations in the eight variables E, H, V, F, 
J, J m , P, and N, which constitute 20 scalars. If the magnetic current density J m is 
specified, equations (1) to (7) form a complete set of equations for the description of the 
dynamics of a warm plasma. 

At this point, the equations will be linearized by expressing the variables as the sum 
of a static term plus a small time -varying perturbation term. Thus, 


E( r, t) = e( r)e^ 

(8a) 

H( r, t) = h(r)e icot 

(8b) 

V( r, t) = v( r)e iu,t 

(8c) 

m (r,t)=J m (7)e lu,t 

(8d) 
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(8e) 


mi i iiiiiii 


N( r, t) = n Q ( r) + n( r)e ia,t 

P( r, t) = p Q ( r) + p( r)e iu,t (8f) 

where a time variation of has been used for all nonstationary quantities. The total 

E, H, V, and J m fields are assumed to consist of only perturbation quantities. The 
quantities n Q ( r) and p Q ( r) represent the equilibrium electron density and pressure, re- 
spectively. 

Substituting equations (8a) to (8f) into (1) to (7) yields the following set of equations: 


icon + v • (n Q v) = 0 

(9a) 

io>mn o v = n Q qe - vp 

(9b) 

rP 0 

p = — - n 
n o 

(9c) 

V X h = n qv + icoe^e 
0^ 0 

(9d) 

V X e = -j m - iu>M 0 h 

(9e) 


The factor e* w *' has been suppressed in the preceding equations. In obtaining equations 
(9a) to (9e), all products of perturbation terms have been neglected in order to make the 
equations linear. 

Equations (2) and (8f) force the equilibrium electron pressure p Q ( r) to be a constant. 
Consequently, the equilibrium electron density n Q ( r) must also be a constant since the 
pressure and density are related by an equation of state. The requirement that the equi- 
librium density be uniform is a direct result of neglecting ionization and recombination. 
The random motion of the electrons due to the nonzero electron temperature will distri- 
bute the electrons uniformly in space unless there is a source or sink of electrons. 

p 

The quantity a = yP 0 /n Q m in equation (9c) can be identified as the square of the 

adiabatic acoustic velocity for the electron gas. If the electron gas obeys the equation of 

o 

state for an ideal gas, a is given by 


a 


2 


rkT e 


m 
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where k is the Boltzmann constant, and T 0 is the electron temperature. 

Since the perturbation electron pressure and density are proportional, equations 
(9a) to (9e) can be reduced to the following set: 


o 

iwp + mn Q a V • v = 0 

(10a) 

iwmn Q v = n Q qe - Vp 

(10b) 

V x h = n Q qv + io>e o e 

(10c) 

V x e = -T m - iwM Q h 

(lOd) 

Equations (10a) to (lOd) are the basic equations for a warm plasma. The standard equa- 
tions that describe a cold plasma can be obtained from the preceding set by setting the 
pressure equal to zero. 

The basic equations can be put into a more tractable form by combining equations 
(10a) to (lOd) to give 

V x v x h - kgh = -ioe 0 e p j m 

(11) 

V 2 p + k 2 p = 0 

(12) 

e=- 1 v X h - ^ vp 

We o e P w2e o e p m 

(13) 

V = - q V X h + 1 vp 

o, 2 £o£p m 

(14) 


2 2 

The new symbols introduced are = n a /cm, the square of the electron plasma fre- 
2 2 poo 

quency, e = 1 - u> p /u> , the effective relative dielectric constant of the plasma, 

,2 2 / 2 , , 2 2 /2 

k e = e p w / c , and k p = e p w /a . 


The previous equations show that the solution for the field in a warm plasma medium 
can be obtained if j is known by first solving equations (11) and (12) for h and p and 
then using equations (13) and (14) to compute e and v. 
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THE MODEL 


The plasma layer surrounding a reentry ve- 
hicle has a very complex structure. The electron 
density and pressure possess both spatial and 
time variations. In order to obtain a model that 
can be analyzed exactly, several simplifying 
assumptions must be made. 

Figure 1 shows the model that will be con- 
sidered. The surface of the space vehicle will be 
approximated by a plane of infinite extent and of 
infinite electrical conductivity. The case where the plane is finite in length has been 
treated by Oliner and Tamir (ref. 8) using the cold plasma model. Their results indicate 
that, as long as the source of the electromagnetic radiation is at least 10 wavelengths 
from the nearest edge of the plane, the infinite plane will represent a valid approximation 
of a finite plane. In practice, the surface of a space vehicle is usually cylindrical. How- 
ever, if the radius of curvature of the surface is much larger than 1 wavelength, the in- 
finite plane model should give useful results. 

The source for the electromagnetic radiation will be a uniformly excited slot antenna 
of width d and of infinite extent located on the conducting plane. The infinite extent of 
the antenna simplifies the problem by making it two dimensional. 

The plasma layer surrounding the reentry vehicle will be approximated by a plasma 
slab of uniform density in the region 0 < x < b. The remaining region x > b is assumed 
to be a vacuum. A plasma slab of this form could never exist in practice unless some 
confining mechanism were provided at the discontinuity in the electron density along the 
surface x = b. However, the slab configuration does provide a useful mathematical 
model. 

In the previous section, it was shown that, in a warm plasma medium, the equilib- 
rium electron density must be constant if ionization and recombination are absent. The 
most general equilibrium density variation that can be used in conjunction with the pre- 
viously developed equations is one that is uniform in space except for possible step dis- 
continuities along a surface. Therefore, the most general plasma slab that could be con- 
sidered would be multilayered. The thickness, density, and temperature of each layer 
could be specified to give the best approximation of the actual density-temperature pro- 
file. However, in the interest of simplicity, only the single-layered plasma slab will be 
considered in this analysis. 

The electrical characteristics of the model can be completely specified by the dimen- 
sionless parameters o^b/c, cod/c, a/c, and In practice, the numerical values 

of these parameters will have wide variations that depend on the configuration of the 
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space vehicle and its reentry trajectory. Typical ranges of values of these parameters 

for space vehicles such as Mercury, Gemini, and Apollo are 0 < cj b/c < 20, 

O P 

0<a/c<2xl0 , and 0 < ci>p/u> < 5. The value cd^b/c = 20 corresponds to a plasma 

layer thickness of 3A n , where A is the free-space wavelength at the plasma frequency, 

p _ 3 p 3 o 

and the value a/c = 2x10 corresponds to an electron temperature of 15x10 K. A 

typical value for u> d/c is 0. 5, which corresponds to a slot width of 0. 1 a q , where a q 

is the free-space wavelength at the operating frequency. 


BOUNDARY CONDITIONS 

Before a formal solution for the field can proceed, it is necessary to select a set of 
boundary conditions for the field components. The boundary conditions that are selected 
should satisfy the physics of the problem as closely as possible and also require the solu- 
tion for the field to be unique. A systematic method for making this selection is first to 
determine all the possible sets of boundary conditions that will yield a unique solution and 
then choose the set that conforms to the physical problem as closely as possible. 

First, consider the plasma region 0 < x < b. The magnetic field in this region sat- 
isfies the differential equation (eq. 11 ), which is repeated here for convenience: 

VXVXh-k^h = - iwe o e p^m (11) 

It is assumed that j m is known. Let f be an arbitrary vector function. Scalar multi- 
plying equation ( 11 ) by f gives 

1 . vx vxh - k^f • H= -iwe 0 e p f • J m (15) 

By using the vector identity 

T • v x v x h = (vxh) • (vxf) + v- [(vxh)xT] 

equation (1 5) can be rewritten in the following form: 

(VXh)- (vxf) - k^F- h + v ■ [(vxh)xf]= icoe^pf * ^m ( 16 ) 

Let both h^ and hg be solutions of equation (11). Then h d = - hg is a solution of 

the associated homogeneous equation. Thus, 
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j|c »j|( 

Now let f = h d (h d = complex conjugate of h d ) and integrate equation (17) over the volume 
of the plasma V : 

f |vxh d | 2 dV-k 2 f |h d | 2 dV + / (vxh d )- (^x^O (18) 

v p p p 

The divergence theorem has been used to reduce one of the volume integrals to a surface 

integral. The vector n is the unit normal to the surface S„, which bounds V . 

2 P P 

At frequencies below the plasma frequency, k g is negative. If the boundary condi- 
tions are selected so that the surface integral is positive definite, it must follow that 
h d = 0 or h^ = h 2 , since equation (18) consists of a sum of positive definite quantities 
equated to zero. Thus, the uniqueness of the magnetic field is ensured provided that the 
surface integral is positive definite. Two obvious choices for making this integral vanish 
would be to specify either nXh or n x (V x h) on the surface of the plasma. 

The boundary conditions for the uniqueness of the pressure field can be derived in an 
analogous manner. Let p d = Pj - P£> where both p^ and p 2 satisfy the differential 
equation for the pressure (eq. 12). By performing the obvious manipulations, it can 
be shown that 



2 

At frequencies below the plasma frequency, k p is negative. The boundary conditions 
that make the surface integral negative definite will require the pressure field to be 
unique. Two possible choices would be to specify either p or 3p/ 3n on the surface of 
the plasma. 

The previous discussion shows that the field will be unique provided that the boundary 
conditions are selected so that the surface integrals that appear in equations (18) and (19) 
are positive and negative definite, respectively. The selection of a set of boundary con- 
ditions from the many possible sets revealed by uniqueness considerations is determined, 
of course, by the physics of the problem. The usual boundary conditions for the electro- 
magnetic field should apply in this problem. Along the surface of the perfect conductor, 
the tangential component of the electric field must vanish It is readily apparent from 
equation (13) that the condition n X e = 0 at x=0 is equivalent to 


( 20 ) 


nx(vxh) 


|x=0 


jq 

com 


n X vp 


|x=0 


where n = a . At the plasma- vacuum interface, the tangential components of both the 
electric and the magnetic field vectors must be continuous. Again, from equation (13), 
the requirement that both n x e and n X h be continuous at x = b is equivalent to 


nX(VXh) - -H- n X vp 
com 


= e n X (v X h) 

r 


x=b" 


x=b 


( 21 ) 


and 


n x h 


= n x h 


x=b" 


x=b"* 


( 22 ) 


where n = a x> Thus, the boundary conditions (20) to (22) will specify n x (v x h) on the 
surface of the plasma region provided that both the pressure field in the plasma and the 
electromagnetic field in the vacuum are known. It has already been shown that specifying 
n X (v x h) on the surface of the plasma is a sufficient condition for the uniqueness of the 
magnetic field within the plasma. 

The boundary conditions for the pressure are more difficult to select since the model 
under consideration is quite idealized. The possible choices of boundary conditions, how- 
ever, can be grouped into two categories. In the first group are those boundary condi- 
tions that require the perturbation pressure to vanish on the surface of the plasma. This 
condition along with the boundary conditions for the electromagnetic field allow the mag- 
netic and pressure fields to uncouple. Since the magnetic current does not appear as a 
source term in the pressure equation (eq. 12), the pressure field will not be excited. 
Therefore, this boundary condition for the pressure will simply yield the solution for the 
cold plasma model. All the other boundary conditions for the pressure will involve cou- 
pling between the magnetic and the pressure fields. 

The boundary condition that will be used is 


/> 

n 




= 0 


lx=0, b” 


(23) 


with n = a x< This boundary condition is obviously an approximation. It has, however, 
been widely used in the literature, which allows this analysis to be compared more easily 
with previous work. 

Equation (23) can be put into the form 
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( 24 ) 


IP 

9n 


in a _ 

n • V X h 


lx=0, b' 


we. 


x=0, b‘ 


by using equation (14). The boundary condition (24) specifies Sp/3n on the surface of 
the plasma provided that the magnetic field is known. It has been shown that specifying 
dp/dn on the surface is a sufficient condition to ensure the uniqueness of the pressure 
field. 

The previous discussion indicates that the boundary conditions that require the sur- 
face integrals in equations (18) and (19) to vanish are sufficient conditions for the unique- 
ness of the field within the plasma. The boundary conditions (20) to (22) and (24) make 
the surface integrals vanish over the surfaces x = 0 and x = b of the plasma slab. How- 
ever, to obtain physically meaningful solutions, the surface bounding the edge of the 
plasma slab at infinity must also be considered. The proper boundary condition to im- 
pose on this surface is a power radiation condition; that is, all waves at infinity should 
transport energy away from the source. A power -radiation condition must be imposed 
in lieu of the usual phase radiation condition since a plasma slab can support backward 
waves. A backward wave has its phase velocity and group velocity vectors in opposite 
directions so that a point on the wave of constant phase may be moving in the negative 
z-direction, for example, whereas the energy carried by the wave is transported in the 
positive z-direction. The power- radiation condition simply ensures that the net energy 
transport at infinity is directed away from the source. 

So far, uniqueness conditions have been derived for the fields in the plasma region 
when the frequency is below the plasma frequency. At frequencies above the plasma fre- 
quency the proof is no longer valid. Resonant modes of arbitrary amplitude may exist 
when the relative dielectric constant of the plasma is positive. However, if an energy- 

loss mechanism such as collisions between the electrons and ions were introduced, the 
2 2 

coefficients k and k in equations (18) and (19) would become complex with negative 
® P — 

imaginary parts for all frequencies. This requires p^ = 0 and h^ = 0, which ensures 

the uniqueness of the magnetic and pressure fields in the plasma. 

The remaining point to consider is the uniqueness of the electromagnetic field in the 
vacuum region x > b. The magnetic field in this region satisfies the equation 

V x v x h - k^h = -ita>e 0 l m (25) 

2 2 2 

where k Q = oj /c . (The equations that the field must satisfy in the vacuum region, 
i. e. , Maxwell's equations, can, of course, be obtained from any of the previous sets of 
equations by systematically setting the equilibrium electron density equal to zero. ) 
Following a procedure identical to that used to obtain equation (18) gives 
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(26) 


L 


VX h. 


dV - k 


2 f 

o 


dV + <f (V x h d ) • (h^ x n)dS = 0 


where V y is the volume of the vacuum region x > b, and S v is the boundary surface, 
which includes the plasma-vacuum interface and the surface at infinity. Since both 
n x h and n x (v x h) are specified on the plasma-vacuum interface by equations (21) 
and (22), this portion of the surface integral will vanish. A sufficient condition to make 
the integral over the surface at infinity vanish is to impose the two-dimensional phase- 
radiation condition 


Urn r 1 / 2 



= 0 


r 2 2~i^/^ 

where r = (x - b) + z for the model under consideration. This condition requires 

the field at great distances from the source to represent a divergent traveling wave. The 

phase-radiation condition is equivalent to a power -radiation condition in this case since 

backward waves no not exist in the vacuum region. If small losses were introduced into 

o 

the medium, the imaginary part of k Q in equation (26) would be negative at all frequen- 
cies. Thus, h d = 0, which ensures the uniqueness of the magnetic field in the vacuum. 

The previous discussion indicates that the following boundary conditions are suffi- 
cient conditions for a unique solution and are those that will be used in the analysis: 

(1) The component of the electric field tangent to the plasma surface must vanish at 
x = 0 and be continuous at x = b. 

(2) The component of the magnetic field tangent to the plasma surface must be con- 
tinuous at x = b. 

(3) The component of the electron fluid velocity normal to the plasma surface must 
vanish at x = 0 and x = b. 

(4) The total field must satisfy the proper phase or power- radiation condition at in- 
finity. 


FORMAL SOLUTION FOR THE FIELD 

The formal solution for the field excited by the slot antenna will be determined by 
using Fourier transform techniques. It will be convenient to represent the slot antenna 
with a magnetic current. The field radiated by a uniformly excited slot antenna is iden- 
tical to the field radiated by a uniform magnetic surface current on the aperture of the 
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antenna with the aperture closed by a perfectly conducting wall. The magnetic current 
density j can be expressed as 


j 


m 




(27) 


where I is the intensity of the total magnetic current. 

An examination of the basic equation governing the magnetic field (eq. 11) reveals » 

that the polarization of the magnetic field will be in the direction of the magnetic current 
j , or in the y-direction for the coordinate system chosen. The other components of the t 

magnetic field are not excited since the boundary conditions are invariant with respect to 
the y- and z- coordinates. Thus, the field components of interest will be the x- and 
z- components of the electric and velocity fields and the y- component of the magnetic field. 

The pressure will be a scalar function of the x- and z- coordinates. All other field com- 
ponents will vanish. 

Let c^(x,/3) and c^(x, /3) be the Fourier transforms with respect to z of h y (x, z) and 
p(x, z), respectively: 


<^(x, 0) = J h y (x, z)e lj3z dz 

/ 'c° 

p(x, z)e 1/3z dz 

oo 


(28a) 


(28b) 


The functions h y (x, z) and p(x, z) can then be recovered by the usual inversion integrals: 


h y (x, z) = -i- j£<^(x, j3)e" 1/3z d/3 

(29a) 

P(x, z) = — /^(x,/3)e-^ z d/3 
2n J C 

(29b) 


where C is a suitable contour in the complex j3-plane. It is assumed that both the 
plasma and vacuum media have small losses so that the fields decay as |z| — °°. This 
is necessary to satisfy the uniqueness conditions and to allow the Fourier transforms to 
be taken. In obtaining numerical results, however, the losses will be neglected. 

The Fourier transform with respect to z of the basic equations for the magnetic 
field and pressure (eqs. (11) to (14)) are easily seen to be 


fr 


I 
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(30a) 



0 < x < b 


3^ 


3x 


(kp - I3 2 )^= 0 


0 < x< b 



b < x 


(30b) 


(30c) 


= 0 b < x 


(30d) 


The right side of equation (30a) is recognized as being proportional to the Fourier trans- 
form of the magnetic current density given by equation (27). The boundary conditions 
given by equations (20) to (22) and (24) transform in a similar manner: 




= o 


3x com 


x=0 


9 qj3 <p 

3x com 



dJf 

3x 


tx=b^ 


(31a) 


(31b) 


& _ n _M\ = 


0X 


n o q ^ 




coe. 


= 0 


x=b" 


(31c) 
(3 Id) 

(31e) 


In addition to the previous boundary conditions, the phase- or power- radiation conditions 
at infinity must also be imposed to specify the solution uniquely. 
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It will be convenient to include the effect of the magnetic surface current in the equa- 
tions for the boundary conditions rather than in the differential equation for the magnetic 
field. This procedure is possible since the source for the field is located along one of the 
boundaries. Consider the surface current to be slightly elevated from the plane x=0 so 
that 6(x) is replaced by 5(x - e) in equation (30a) where e « b. Since the magnetic 
field is continuous across the magnetic surface current, e^(e", /3) is equal to <^(e + , $). 
Therefore, the boundary condition (3 Id) is unchanged. The derivative of the magnetic 
field with respect to x is discontinuous, however, at x = e by the amount 


^oVm 



Thus, equation (31a) can be modified to read 


dcJf _ q/3 tp 
dx com 


sin 


/3d 


^oVm 


2 


(32) 


The solutions for the transformed variables Jl and ^ can now be determined where 
equation (30a) is regarded as homogeneous. 

A general solution of equations (30a) to (30d) is given by 


/ -ilYx 
Hje 


+ R^e 


ir h x ^ 


0 < x < b 


(33a) 


£P = P 


/ -ir x ir x 
e p p 


0 < x < b 


(33b) 


Jf=H Q (e irht> + R 1 e irhb je" Q '^ x " b ^ b < 


= 0 b < x 


(33 c) 
(33d) 


where 



(34a) 



(34b) 
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a = [ p 2 - k 2 


A 1/2 


o/ 


(34c) 


and H q , P q , Rj, and Rg are complex functions of /3 to be determined from the bound- 
ary conditions. The boundary condition that requires the magnetic field to be continuous 
at x = b has been built into the general solutions. Substituting the general solutions into 
the remaining equations for the boundary conditions yields the following system of equa- 
tions: 


-i/3 

-i/3 

o>€_r_ 

o p 

n 0 q 

-coe r 
op 

n 0 q 


H o 


0 

-ilYb 
-ij3e h 

ir.b 
-i/3e h 

-ir b 

o P „ P 

Vi 

- W£ o r p e ir p b 

V* 


H o R l 


0 

ir h 

- ir h 

q£ 

sL 


p 


sinM 
-ico^ * I ^ 

com 

com 


0 


o pm ^ 








2 

-ilYb 

< ir h - “V e 

ir.b 

-( ir h + “ e p )e 

5Le' ir P b 

com 

-3£ e ir P b 

com 


P o R 2 


0 


The solution for the transformed field components <A and cP can be determined by 
solving the preceding system of equations for H Q , H Q R^, P Q , and P Q R 2 , and then sub- 
stituting the results into the general solutions (33a) to (33d). The result of this operation 
is 


D(/3) 

0 < x < b 

(35a) 

S(x,p) - M ( X ^) 
D(/3) 

0 < x < b 

(35b) 

J(x,0) = NM e-“ 
D(J8) 

(x- b ) b < x 

(35c) 

#(x,/3) = 0 

b < x 

(35d) 


where 
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/3d 


N(x,0) = 


,3 -p j sin 
e o e P r P I m 2 


0) 


£d 

2 


I/3 2 sin r h x + /3 2 sin T h (b - x)cos r b 


2 2 

+ — r h r p c os r h (b - x)sin r p b + ff£ p — T p sin T h (b - x) sin T p b 
2 2 


co. 


CO. 


(35e) 


M(x,/3) = 


ia,2e o e P m ^m 


sin 0d 
2 


2 w 

/3 sin r\b sin r (b - x) r IY cos r x 

n P 2 P n P 


co. 


2 2 
— r r. cos IYb cos r (b - x) - ae — T sin r.b cos r (b - x) 

2 P n n p PoP 11 P 


co. 


P2P 

co 

P 


(35f) 


DO) = 2/3 2 r h r p ^ (l - cos r h b cos r p b) + (b 4 + r 2 r 2 i^Wn r h b sin r p b 

W D \ “V 


ae p r P ~2 r Sin r h b C0S r p b + r P Fh \ 1 C0S r h b Sin r P b (35g) 

co p \ 


The solution for the transformed field is not unique until the P ro P er branches of the 
double valued P arameters T h , T p , and a, as defined by equations (34a) to (34c) are 
specified. Since both and & are even functions of and T p , the choice of either 
branch for these two parameters will yield the same result for ^ and The branch 
of a where Re a ^ 0 and Im a ^ 0 must be selected to ensure that Jl either van- 
ishes at infinity or corresponds to an outward propagating wave. This requirement on a 
is a consequence of the phase -radiation condition. 

The solution for h y and p can now be determined by employing the inversion inte- 
grals (29a) and (29b) where the contour C is shown in figure 2. For z > 0, the contour 
can be closed with a semicircle of infinite radius in the lower half of the /3-plane. The 
contour must be deformed around the branch cut to remain on the proper branch of a. 
The contour must also be deformed around the poles of the integrands along the real 
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Re P 


Figure 2 . - Complex f^plane. 


/3-axis. The only poles that should be encir- 
led by the closed contour are the poles that 
correspond to waves that transport energy 
away from the source at infinity. This re- 
striction on the contour is imposed by the 
power- radiation condition. Thus, the forward 
wave poles located along the positive real 
/3-axis and the backward wave poles located 
along the negative real /3-axis should be en- 
circled. The nature of the poles (i. e. , 
whether they correspond to forward or back- 
ward waves) can be determined by computing 
the phase velocity 
velocity v„ 


g 


Vp = oj// 3 and the group 
= 9w/9/3 of the associated waves. 


If the product VpV^ is positive, the wave is 
a forward wave; whereas, if it is negative, the wave is a backward wave. 

It can be shown that the integration over the portion of the contour C M vanishes. 
Therefore, from Cauchy's theorem, the integral along the contour C is equal to -2id 
times the sum of the residues of the enclosed poles plus the branch cut integral. Thus, 


h y (x,z) = 



iN(x,/3 n ) 


9D(/3) 

9/3 


-i/3 z 

e + Branch cut integral 0 < x < b (36a) 


/3=/3 n 



iM(x,/3 ) -i/3 z 

p(x,z) = ) . e + Branch cut integral 0 < x < b (36b) 


n=l 


3DQ) 

9/3 


/3=/3 n 


hy(x, z) = 



-a (x-b) -i/3 z 

5 e + Branch cut integral b < x (36c) 


21 


A 

I! 



p(x,z) 3 0 b < x 


(36d) 


where the j3 n are the roots of D(/3) = 0. The terms in the summations in equations 
(36a) to (36c) correspond to the discrete spectrum of waves guided by the plasma layer 
that are often referred to as surface waves, whereas the branch cut integrals correspond 
to the continuous spectrum of waves of the radiation field. Note that the discrete spec- 
tra of h y and p are identical, as they must be, since the magnetic and pressure 
fields are coupled by the boundary conditions. 

A similar procedure can be carried out to evaluate h^ and p for z < 0. For this 
case, the contour would be closed in the upper half of the /3- plane. This completes the 
formal solution for the field. In order to proceed with the analysis it will be necessary 
to solve for the roots of D(/3) = 0 and to evaluate the branch cut integrals. These topics 
will be covered in the following two sections. 


DISPERSION CURVES 

The roots of the equation D(/3) = 0 can be divided into two groups. In the first group 
are the roots that correspond to the choice Re a > 0. These roots generate physically 
realizable solutions for the field; that is, they satisfy the radiation condition. The sec- 
ond group contains the roots corresponding to the choice Re a < 0. These roots generate 
solutions that are not physically realizable since they ’’blow up” at infinity. In this anal- 
ysis, only the roots in the first group are of interest. 

The transcendental equation D(/3) = 0 has both real and complex roots. The loci of 
these roots in the complex j3-plane possess a great deal of symmetry. If (3 = B is a 
root, /3 = -j3 q is also a root since the variable ft always appears in D(/3) as /3 . This 
phenomenon should be expected since the plasma medium is isotropic. If a solution is 
obtained that corresponds to a wave traveling in the positive z-direction, there should 
also be a similar solution that corresponds to a wave traveling in the negative z-direction. 
In addition, if = /3 Q is a root of D(/3) = 0, /3 = /3* is also a root. Because of this sym- 
metry, it is only necessary to determine the roots of D(/3) = 0 in one quadrant of the 
complex /3-plane. 

If one is interested in the energy transport along the plasma layer, it is only neces- 
sary to compute the real roots of D(/3) = 0. Some of the complex roots may correspond 
to physically realizable solutions. However, the waves associated with these complex 
roots cannot transport any net energy. These roots would only be of interest if a calcula- 
tion of the field pattern in the immediate vicinity of the slot antenna were desired. It 
should be noted that, if the original equations for a warm plasma contained dissipative ef- 
fects such as electron- ion collision terms, all the roots of D(/3) = 0 would be complex. 
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In order to eliminate the tedious numerical work involved in computing the complex zeros 
of the transcendental function D(/3), dissipative effects were neglected. 

The real roots of D(/3) = 0 were computed on an IBM 7094 computer for various fre- 
quencies. It was found that at frequencies below the plasma frequency there are, at most, 


four real zeros of D(/3). However, at frequencies above the plasma frequency, the num- 

/ 2 / 2 ^ ^ 

ber of real zeros is usually quite large, of the order of (cob/ ira.) (1 - ajp/w J . To avo: 
dealing with such a large number of zeros and still cover a reasonable frequency range, 


To avoid 


unrealistic values must be used for either the plasma layer thickness or electron temper- 


ature. Qualitative results for more realistic values of these parameters can then be 


easily inferred. 

2 2/2 

The computations were performed for the parameters a) b / c =5 and a/c = 0. 1. 
22/2 P ' 

The value of a)pb /c corresponds to a plasma layer thickness of 0. 35 Ap, where Ap 

is the free- space wavelength at the plasma frequency, and the value of a/c corresponds 

to an electron temperature of 37X10^ °K. The real roots of D(/3) = 0 are shown in fig- 

ty 

ure 3(a) in the form of dispersion curves with as a function of (co/u>^) , where 



(0/ Up ) 2 


(b) Cold plasma. 

Figure 3. - Dispersion curves for wave coupled to grounded plasma layer. 
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X Q is the free-space wavelength of the waves, and X g is the wavelength of the guided 
waves in the plasma layer. (The curves have been numbered for convenience.) Since all 
the curves are located in the range 1 ^ X^Xg < c/a, it can be concluded that the phase 
velocities associated with the waves are bounded from above by the velocity of light and 
from below by the electron acoustic velocity. Curves 3 to 8 in figure 3(a) appear to start 
at points along the line \,Ag = 1* However, it must be remembered that these curves 
show only the loci of the real roots of D(/3) = 0 with frequency. At frequencies below 
the plasma frequency, these roots are complex, which causes the associated waves to be 
cut off. As the frequency increases, these roots will move in the complex /3-plane to- 
ward the real /3-axis, and at critical frequencies for each root, they will reach the real 
/3-axis at the point /3 = k Q . As the frequency further increases, they will rapidly move 
down the real /3-axis away from the origin and asymptotically approach the point /3 = co/a. 

The dispersion curve for the corresponding cold plasma case is shown in figure 3(b), 

where the cold plasma layer is shown to support only one propagating wave and then only 

o 

in the frequency range (co/co p) < 0. 5. This result is in sharp contrast with the warm 
plasma layer that can support several propagating waves at all frequencies. A compari- 
son of figures 3(a) and (b) reveals that curve 2 in the warm plasma case corresponds to the 
single curve in the cold plasma case in the frequency range (co/co ) < 0. 5. The remain- 
ing portion of curve 2 and all the other curves in figure 3(a) are, therefore, unique to the 
warm plasma layer. 

The transverse variations of the magnetic and the pressure fields for the different 
guided waves (fig. 4) can be determined from the expressions for N(x,/3) and M(x,/3), 
respectively. It can be shown that the wave associated with curve 1 in figure 3(a) is very 
tightly bound to the plasma- conductor interface; that is, both N(x,/3j) and M(x,/3j) de- 
crease rapidly in amplitude with increasing 
x. This wave is essentially the same wave 
as derived by Seshadri (ref. 9) for the case 
of a semi- infinite warm plasma bounded by 
a conducting surface. The wave that cor- 
responds to curve 2 is tightly coupled to the 
plasma- vacuum interface. All the field 
components associated with this wave decay 
rapidly in amplitude as jx - bj deviates 
from zero. This wave is essentially the 
same wave as derived by Hessel, et al. 

(ref. 10) for the case of a semi- infinite 
warm plasma bounded by a vacuum. The 
waves associated with curves 3 to 8 are a 
result of an interaction between the plasma- 



wave l 

Plasma layer 
thickness, b 



Rectangular coordinate, x — 

Figure 4. - Typical transverse variations of waves. 
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conductor and plasma- vacuum interfaces. The magnetic field intensity associated with 
these waves has its maximum value at either x = 0, the plasma- conductor interface, or 
x = b, the plasma-vacuum interface, and it is quite small between these surfaces, where- 
as the pressure has a standing wave pattern between these two surfaces. 


FAR-ZONE RADIATION FIELD 

The branch cut integrals give rise to the radiation field. An asymptotic evaluation 
of the far- zone field will be determined by using the method of saddle-point integration 
(ref. 11). The solution for the transformed magnetic field in the region x > b is given 
by equation (35c), which is repeated here for convenience: 

(35c) 

om 


The magnetic field 


is then given by 


h(x,z)=-L f e ' a ^e 1/32 d/3 x>b (37) 

y 2 ttJ q D(/3) 


where the inversion contour C is shown in figure 2 (p. 21). 

At this point, it is convenient to map the complex /3-plane into the complex ^-plane 
with the mapping function 


v 



/3 = k Q sin v (38) 

where v = a + itj, and to transform 
x and z into their cylindrical coordi- 
nate equivalents with the transforma- 
tion z = r cos 6 , x = b + r sin 0 (see 
fig. 1, p. 10). The mapping of the 
/3-plane into the y-plane is shown in 
figure 5. The entire /3-plane maps 
into a strip of width 7 t in the y-plane. 
The quantities Ql, Q2, Q3, and Q4 
refer to the images of the quadrants of 
the /3-plane in the i'-plane. With these 
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transformations, the expression for the magnetic field becomes 




N(b,k Q sin i^cos v e 


-ik Q r sin(7>+0) 


dv 


D(k Q sin v) 


( 39 ) 


The integral for h y (r, d) will be evaluated by the method of saddle-point integration. 
This method consists of deforming the contour C into a new contour C* so that the 
magnitude of the exponential term in equation (39) is very small everywhere except in the 
vicinity of one point along C\ The exponential term will behave much like a delta func- 
tion in that the value of the integral will be proportional to the value of the integrand at 
this one point. 

The saddle point of the exponent -ik o r sin(y + 9) occurs where cos(v + 9) = 0 or 
where v = 77/2 - 9. In the vicinity of the saddle point v = v/2 - 9, the function 
-ik Q r sin(^ + 9) can be expanded in a Taylor series to give 

k r « 

-ik Q r sin(^ + 9) = -ik Q r pr sin <p cos cp + . . . 


where pe 1 ^ * v - (n/2 - 9). The real part of -ik Q r sin(y + 9) will have the greatest rate 
of change along the contour that passes through the saddle point at an angle of cp = 77 / 4 . 
This steepest descent contour (SDC) is shown in figure 6. 

The exponential term in equation (39) has a unit modulus at v = n/2 - 9 and falls off 

9 

as exp[-(k Q r/2)p ] with distance p along the steepest descent contour. If k Q r » 1, 
the dominant portion of the integral will be the integrand evaluated at v = v/2 - 9 times 
the integral of the Gaussian term with respect to p. Thus, 


k N(b,k cos 0)sin 9 -i(k r-7r/4) 

h v (r, 9) e 

y 277 D (k Q COS 9) 


/ 

‘'SDC 


(-k r/2)p 

e 0 dp (40) 


In evaluating the integral in equation (40), the limits of integration can be selected as 
-00 and +00 since the major contribution to the integral comes from the portion of the con- 
tour near p = 0. Performing the integration gives 


hy(r, 9) 


k o N < b > k o 
2tt 


-i(k r- 77 / 4 ) 

cos 0)sin 9 e 


D(k Q cos 9) 



1/2 


(41) 
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If more terms were retained in the 
expansion of -ik Q r sin(^ + 0), the ex- 
pression for h . would contain addi- 
y 3/2 

tional terms of the order (k Q r) _ / and 

lower. The additional terms would give 
a better approximation to the magnetic 
field hy especially if k Q r is not ex- 
tremely large. However, these addi- 
tional terms do not contribute to any net 
radiated power. 

In deforming the contour C into 
the steepest descent contour, some of 
the poles of the integrand may be 
crossed. This usually happens if the 
observation angle 6 is near 0 or ir. 
interface, the total field also consists 

of the discrete spectrum of waves. 

An interesting phenomenon occurs if one of the poles corresponds to a backward 
wave. In this case, the steepest descent contour encircles the backward wave pole twice, 
as shown in figure 6. If k Q r is very large, the effect of these backward wave poles will 
not be noticed. However, if k Q r is small, the location of the zeros of D(k Q sin v) in 
equation (39) will have an appreciable effect on the field. The backward wave poles will 
give rise to a standing wave since two poles are encircled. This standing wave will only 
be observed for observation angles in the vicinity of tt/2. As the observation angle de- 
creases, a point will be reached where only one of the backward wave poles will be en- 
circled. This is the beginning of the region where the backward waves correspond to 
traveling waves. 

Finally, it should be noted that the steepest descent contour crosses the branch cut 

/ 2 2 

in quadrants 2 or 4 of the /3-plane. On the improper branch of 1/3 - k Q I , the contour 
may encounter the improper or nonspectral roots of D(k Q sin v) = 0. The improper roots 
usually give rise to the dominant portion of the near- zone field (ref. 5). These roots 
could be used, for example, to evaluate the field along the plasma- vacuum interface. The 
far- zone field could then be calculated by using a Kir chhoff- Huygens integration over the 
plasma- vacuum interface. However, the extreme difficulty in obtaining the complex 
roots of D(k Q sin v) = 0 precludes the use of this method. 

The shape of the radiation pattern of the far- zone field will be dependent on both the 
parameters that describe the warm plasma layer and on the width of the slot antenna. 
Since the modification of the radiation pattern by the warm plasma layer is of interest, it 
will be expedient to set the width of the slot antenna equal to zero. A slot antenna of zero 



It simply indicates that, near the plasma- vacuum 
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width, which is nothing more than a line source, radiates an isotropic field in the absence 
of a plasma layer. Thus, any deviation from an isotropic pattern when the warm plasma 
layer is present is entirely a result of the plasma. 

The general features of the radiation pattern of the far- zone field can be determined 
from a study of the functions N(b,k Q cos 9) and D(k Q cos 9) appearing in equation (41). 

It can be shown that both N(b,k Q cos 9) and D(k Q cos 9) are strongly dependent on 
and Tp where 



2 2 \ 1/2 
- k“ cos^ 9) 


and 


r 

P 


- *o cos2 S ) 


1/2 


In particular, the radiation is intense for those directions in space where and r 

are real, and it is weak for those directions where and Tp are imaginary. 

At frequencies below the plasma frequency k^ and k^ are negative. Since both 

e P 

T h and Tp are imaginary for all angles 9 in this frequency range, the radiation will 

be relatively weak. At frequencies above the plasma frequency is real for the angles 

2 2 2 2 2 2 
where k - k cos 9 > 0 and IT is real for the angles where k - k^ cos 9 > 0. 

e o Pi 1/ \l/2 P 0 

Thus, there is a critical angle # c = cos" k e /k Q = cos” ( € D ) c/a ass °ciated with 

and a second critical angle 9 Q = cos - ^kp/k Q = cos”*^^*/^ associated with Tp. The 

radiation will generally be weak where j cos 9 J > | cos 9 Q | , | cos 9 1 > | cos 9 Q | and 
strong where |cos 0| < |cos0 c |, |cos0|< | cos 0 |. 

At frequencies above the plasma frequency, the critical angle 9 Q is always located 
in physical space since 0 <( € p)^^ ^ 1* The critical angle 0 Q , however, is usually 
located in some invisible portion of space since (e c / a is usually greater than 1 un- 
less the frequency is very close to the plasma frequency. Since (e ^ c/a is usually 
much greater than 1, variations in the acoustic velocity a will not transform T p from 
a real to an imaginary quantity. Therefore, variations in the electron temperature 
should not significantly alter the radiation pattern unless the source frequency is very 
close to the plasma frequency or the electron temperature is extremely high so that 
a/c ~ 1. These limiting cases have been treated by Stewart and Caron (ref. 6). 

Some typical radiation patterns are shown in figure 7 for the frequencies 
(w/cOp^ = 0.1 and 2. 0. The radiation patterns have been normalized by setting the in- 
tensity in the direction 9 = tt/2 equal to 1. Because of the symmetry of the structure, 
only the range of angles 0 < 9 < 7t/ 2 has to be considered. 

The radiation pattern is not affected by variations in the electron temperature as 

o 

long as the temperature is sufficiently low so that a/c ^ 10 except for the special case 
where the source frequency is extremely close to the plasma frequency. For higher tem- 
peratures, the pattern is slightly changed, as shown in figure 7. The shapes of the 
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Figure 7. - Radiation pattern of far-zone field. 

radiation patterns can be easily explained from ray theory by using the cold plasma model 
(e. g. , see ref. 3). 

If the width of the slot antenna were nonzero, the intensity of the radiation field in 
the direction 6 would be reduced by the factor 


/k d N 
sinl cos 0 


k d 


cos 0 


The reduction in field intensity in all directions is not large unless k Q d exceeds 1. A 
typical value for k Q d/2 is of the order 0.2, which corresponds to a slot width of 0. 1 
X Q , where is the free- space wavelength at the source frequency. Thus, it can be 
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concluded that the radiation patterns, as shown in figure 7, are typical of patterns that 
are associated with a plasma- covered slot antenna. 


ENERGY TRANSPORTED BY THE FIELD 


The final topic to be considered is the energy transported by the field. In an ordinary 

dielectric medium, the energy transport can be computed by using the Poynting vector. 

To determine how the Poynting vector should be modified to include the effects present in 



a warm plasma, consider the divergence of e X h + pv*: 


• (e x h* + pv*) = (v x e) • h* - (vx h*) • 


e + pV 


_* 

v 


(r-v); 


(42) 


The right side of equation (42) can be simplified by using the basic equations (11) to (14). 
The result of this simplification is 


— Jf: % 

(e X h + pv ) = -4ico 



e o e 


— * 
e 



(43) 


The terms in the first parentheses on the right in equation (43) are energy storage terms. 
They are the magnetic, electric, electron potential, and electron kinetic energy densities, 
respectively. The last term in equation (43) represents the energy loss due to the mag- 
netic currents. Equation (43) can be put into the following form by using the divergence 
theorem: 


-Re (f 

2 Js 

Since the right side of equation (44) represents the rate at which energy is being dissi- 
pated in the volume V, the left side of equation (44) must be the rate at which energy is 

— J * * 

being supplied to V. Thus, it is proper to interpret the quantity S =-g Re(eXh +pv ) 
as the energy flux vector for a warm plasma medium. The energy flux vector S is 
simply the sum of the usual Poynting vector e X h , which accounts for the electro- 

jjc 

magnetic energy transport, and the term pv , which accounts for the energy transport 
due to the electron motion. 

Next, consider the power orthogonality properties of the guided waves. Assume that 


(e X h* + pv*) • dS = 


-JJ 


h* • ] m dv 


(44) 
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there are two guided waves having the components hp Pp Vp /3j| and 
| eg, hg, Pg, Vg, /3gl that propagate along the plasma layer. The total power associated 
with these waves per unit length in the y-direction is given by 

P total = P 1 + P 2 + P 12 (45a) 


where 


p i = -g Re f (ej X h* + p.v*) • a z dx (45b) 

•'O ' 

and 

P ij =i Re f (®i * + ® j X \* + P^* + Pf3 • dx (45c) 

J Q 

th 

It is obvious that is the total power associated with the i wave if it alone were 
present. It is expected that P^ = 0 if # 0j. The fact that P^ is zero if |3j * j3j 
can be shown in several ways. First, general expressions for e and v in terms of 
h and p could be substituted into equation (45c). The integration could then be carried 
out by parts giving quantities that are evaluated only at the end points. Employing the 
boundary conditions will require these end-point quantities to vanish. 

A second approach to the problem is to recognize that, for the case where there are 
no energy losses due to dissipative effects such as collisions, the total power associated 
with the waves must be independent of the z-coordinate. Since both Pj and Pg are in- 
dependent of z, it must follow that P^g is also independent of z. In forming P^g, 
-i(/3 r 0g)z i(/3 r /3 ) z 

terms like e and e appear in the integrand (for the assumption that 

j3j and /3g are real). Therefore, if P^g is to be independent of z, either the coeffi- 
cients of the preceding terms must vanish or = /3g . However, since /3j * /3g the co- 
efficients of these terms vanish making P 12 = °‘ 

This proof can easily be extended to show that the power associated with the continu- 
ous portion of the spectrum is also orthogonal to the powers of all the discrete waves. 

The continuous spectrum can simply be regarded as the superposition of an infinite num- 
ber of discrete waves. The powers of each of these waves must be orthogonal to every 
member of the discrete spectrum since the range of j3 for the continuous spectrum does 
not include members of the discrete spectrum. Therefore, the total energy transported 
by the field is simply the sum of the energies transported by each member of the discrete 
spectrum plus the energy transported by the waves of the continuous spectrum. 


Curve 



Figure 8. - Relative powers associated with waves. 


The energy flux density vector S was evaluated for each of the guided waves and the 
radiation field. The energy transported by a guided wave was computed by integrating 
the corresponding vector S over a surface normal to the z-axis and extending over the 
range 0 ^ x < <» and 0 s y < 1. This result was then doubled, since for every guided 
wave propagating in the positive z-direction there will be an identical wave propagating 
in the negative z-direction. The energy transported by the radiation field was computed 
by integrating the corresponding vector S over a segment of a cylindrical surface of 
infinite radius extending over the range 0^9 ^ it/ 2 and 0 ^ y ^ 1. This result was 
then doubled, since the pattern of the radiation field is symmetrical about 9 = it/ 2. 

Numerical results for the energy transport are presented in figure 8 for the param- 
eters o^b /c = 5, a/c = 0. 1, and d = 0. It is instructive to consider first the case 
where the slot antenna has zero width (i.e. , d = 0) and later extend the results to include 
the effect of a varying slot width. The power associated with each wave is expressed as 
a ratio in decibels of the wave power to the total power supplied by the source. The 
curves in figure 8 are numbered for easy comparison with the dispersion curves in fig- 
ure 3(a)(p. 23). 

The results show that, at frequencies below the plasma frequency, the dominant 
mechanism for transporting energy is wave 1, the wave coupled to the plasma- conductor 
interface. The radiation field and wave 2, the wave coupled to the plasma-vacuum inter- 
face, each account for about 1 percent of the total power. 

The nature of these results should be expected. The wave coupled to the plasma- 
vacuum interface decays very rapidly with distance from the interface. Consequently, 
there will be very little coupling between this wave and the source. The wave coupled to 
the plasma- conductor interface, however, is very intense at the location of the source 
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and, therefore, is strongly excited. The radiation field is weakly excited since the plasma 
layer is quite thick and opaque to electromagnetic waves in this frequency range. 

At frequencies above the plasma frequency, the radiation field is more strongly ex- 
cited. However, wave 1 still accounts for the major portion of the power until ^u>/ w p 2 
exceeds about 1.7. At the onset of propagation of waves 3 to 8, the relative intensity of 
the radiation field takes a sharp drop. As the frequency increases, the fraction of the 
total energy that is transported by these additional waves drops quite rapidly allowing the 
relative intensity of the radiation field to increase. Even though the number of these ad- 
ditional waves increases quite rapidly with frequency, the relative amount of energy that 
they transport is decreasing. Thus, in the limit of infinite frequency all the energy will 
be transported by the radiation field regardless of the electron temperature. 

In the cold plasma model, the radiation field transports all the energy at frequencies 
where (oj/co^j^ > 0. 5, since no propagating guided waves can exist in this frequency 
range. This is in sharp contrast to the results for the warm plasma model. The radia- 
tion field in the warm plasma model accounts for only 1 percent of the total power sup- 
plied by the source at (w/o> ) 2 = 0. 5 and about 50 percent of the total power at 

(u>/u> p ) 2 = 2. 

It must be remembered that the numerical results are for an extremely high electron 
temperature. However, the results demonstrate that the electron temperature can have 
an appreciable effect on the intensity of the radiation field. A limited amount of numeri- 
cal work has been done by using more realistic values for the electron temperature. For 
the case where u> b /c =4.0X10” and a/c = 2x10” , which corresponds to a plasma 

r 3 O 

layer thickness of 0. 01 and an electron temperature 15X10 K, it was found that 

wave 1, the wave coupled to the plasma- conductor interface, still accounts for the major 

portion of the power supplied by the source in the frequency range ( oj/o )^ 2 < 1 even 

though the plasma layer is extremely thin. Thus, it can be concluded that at frequencies 

below the plasma frequency the wave coupled to the plasma- conductor interface is the 

dominant energy- transport mechanism for the case where the slot antenna has zero width. 

th 

If the slot antenna has a finite width the power of the n n guided wave will be reduced 
by the factor 


sin 


^n d 


B d 
H n 


If the guided wave amplitudes are to be significantly reduced, the slot width d must be 
sufficiently large so that /3 n d/2 = 7rd/A g > it. Since \,A g is of the order c/a for each 
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of the guided waves except in the vicinity of their low frequency cutoff points (see fig. 3, 
p. 23), the condition for a significant reduction in guided wave amplitudes can be stated 
as d/A Q > a/c. Since a/c is usually in the range 10"^ to 10”^, the condition d/A Q > a/c 
is readily satisfied. It has been shown in a previous section that the finite width of the 
source will not distort or attenuate the far- zone field if k Q d < 1. Therefore, it can be 
concluded that for the case of a slot antenna of finite width the dominant energy- transport 
mechanism is the radiation field except at frequencies near the guided wave cutoff fre- 
quencies. 


CONCLUSIONS 

The radiation from a slot antenna in a ground plane covered by a warm plasma layer 
has been determined. The radiation pattern of the far- zone field was not affected by the 
plasma temperature unless the source frequency was near the plasma frequency or the 
plasma temperature was extremely high so that the electrons have relativistic thermal ve- 
locities. The complete spectrum of waves guided by the warm plasma layer was also 
determined. The spectrum of waves associated with a warm plasma layer was signifi- 
cantly different from the spectrum associated with a cold plasma layer. The energy 
transported by the guided waves often greatly exceeded the energy transported by the 
radiation field for the case where the slot antenna had zero width. However, the radia- 
tion field proved to be the dominant energy-transport mechanism for the case where the 
slot antenna had a finite width. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 28, 1966. 
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